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Abstract 

We present a new formulation of the correlated electron-ion dynamics 
(CEID) by using equations of motion for nonequilibrium Green's func- 
tions, which generalizes CEID to a general nonequilibrium statistical en- 
semble that allows for a variable total number of electrons. We make a rig- 
orous connection between CEID and diagrammatic perturbation theory, 
which furthermore allows the key approximations in CEID to be quan- 
tified in diagrammatic terms, and, in principle, improved. We compare 
analytically the limiting behavior of CEID and the self-consistent Born ap- 
proximation (SCBA) for a general dynamical nonequilibrium state. This 
comparison shows that CEID and SCBA coincide in the weak electron- 
phonon coupling limit, while they differ in the large ionic mass limit where 
we can readily quantify their difference. In particular, we illustrate the 
relation between CEID and SCBA by perturbation theory at the fourth- 
order in the coupling strength. 



1 



1 Introduction 



One of the most fundamental problems in molecular electronics is to under- 
stand the inelastic scattering effects of atomic vibrations on transmitted elec- 
trons. These effects have been extensively studied both experimentally and 
theoretically in recent years [IJ. The interplay between electronic and nuclear 
dynamics in atomic-scale devices influences not only the device characteristics, 
e.g. inelastic current-voltage spectroscopy [2J, but also the device stability due 
to local heating within the junction [3]. 

NonequiHbrium Green's function (NEGF) theory |4l|5l|6j provides a system- 
atic framework for describing the effects of the coupling between transmitted 
electrons and atomic vibrations [71 [3 [U [10] . For weak electron-phonon coupling, 
a well-known approximation for evaluating the dressed Green's function is the 
self-consistent Born approximation (SCBA) [HI O [T0| which sums only non- 
crossing diagrams in the diagrammatic perturbation expansion of the Green's 
function. Because the percentage of noncrossing diagrams decreases quickly 
with increasing order in the coupling strength, SCBA breaks down at strong 
electron-phonon coupling [llj . 

Alternatively, correlated electron-ion dynamics (CEID) [2] [HI [13] has been 
developed for describing the effects of the electron-ion correlation and interac- 
tion on the inelastic dynamics of the electrons and nuclei. CEID, as an extension 
of molecular dynamics, reinstates the electron-ion correlation and the quantum 
nature of nuclei in order to take account of energy exchange between electrons 
and nuclei reasonably. So far CEID has been applied to a wide range of trans- 
port properties of atomic wires, including inelastic current -voltage spectroscopy 
|2] , the calculation of local heating (and its signature on the current) in real time 
when combined with electronic open boundaries [H], and the non-conservative 
nature of current-induced forces [H]. Recently, a comparison between CEID 
and the NEGF method in SCBA has been made both numerically and analyti- 
cally for steady state transport [TB]. However, there are two restrictions on the 
original CEID discussed above. First, it is assumed that the electron-ion sys- 
tem is described in terms of an ensemble with a fixed total number of electrons. 
Second, the original CEID methodology lacks a systematic scheme to improve 
its accuracy. 

In this paper, we develop a new formulation of CEID by using equations of 
motion for a set of nonequilibrium Green's functions which are closely linked to 
the dynamical variables in the CEID method. To illustrate this idea, we consider 
a model system of noninteracting electrons linearly coupled to a quantum oscil- 
lator. The motivation behind this effort is to lift the restrictions on the original 
CEID. We attempt to make a rigorous connection between CEID and diagram- 
matic perturbation theory, so as to quantify the key CEID approximations in 
diagrammatic terms, and, in principle, to be able to improve them. Moreover, 
in the framework of NEGF, the scope of CEID can be readily extended to a 
general nonequilibrium ensemble with a variable total number of electrons. We 
compare analytically CEID with SCBA for a general nonequilibrium state in 
the time domain, thus extending the previous comparison p6] for a steady state 
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in the energy domain. 

The paper is organized as fohows. In the next section, we present the new 
formulation of CEID for the model system by using equations of motion for 
a set of nonequihbrium Green's functions. In section [Sj CEID is analytically 
compared with SCBA for the model system at the fourth-order in the coupling 
strength and in two specific limits: weak electron-phonon coupling limit and 
large ionic mass limit. To evaluate the Green's functions for the system of 
electrons coupled to a classical oscillator, a classical version of Wick's theorem 
is introduced in the appendix. Finally, conclusions are drawn in section IH 



2 Model and formulation 

We consider an infinite open system of noninteracting electrons linearly cou- 
pled to a single quantum oscillator. The electrons are described in terms of 
the second-quantized field operators ^'(r) and ^'"'"(r). The Hamiltonian of the 
system then takes the form 



H = y dr*+(r) (^-^ ) *W + 

drF{r)'^+{r)'^{r) 



2M 
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-X 



(1) 



with X = R — Rq. The first two terms constitute the free-particle Hamiltonian 
Ho, and the last term describes the electron-phonon interaction W. Here V{r) 
is the lattice potential and F{r) is the electron-phonon couphng strength. Rq 
and K are the equilibrium position and the spring constant of the harmonic 
oscillator respectively. 

In the absence of electron-phonon interaction if* (taken to exist at t = 
— oo), the unperturbed electron subsystem was settled in the Landauer steady 
state, which is characterized by two sets of one-electron states, i.e. Lippmann- 
Schwinger scattering states, {|$iQ)} (a = 1,2) with occupancies fia set by the 
battery terminals [iTt I18j. In the |$iQ)-representation, the statistical operator 
of the unperturbed system is thus taken to be 



Po = ^ exp ■ 



\2M 



-KX' 



where (3 = l/ksT is the inverse temperature, electrons occupying the two sets 
of Lippmann-Schwinger scattering states are characterized by the chemical po- 
tentials fia {a = 1,2) respectively, and Z is a normahzation factor ensuring that 
Tr(po) = 1- Here, {c^} and {cia} are the creation and annihilation operators 
for the complete and orthonormal set of the Lippmann-Schwinger scattering 
states The fermion field operator ^'(r) (^'+(r)) can thus be expressed 

as a linear combination of {cia} ({c^})- 

We now define the contour-ordered Green's function 
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Girt, r't') = [ih)-^ {Tci^HirmWt')) (2) 

where the angular bracket (• • •) = tr (po ■ • ■)• By virtue of the grand-canonical 
structure of po, we have thus allowed for an ensemble with a variable total 
number of electrons. The contour C runs from t = — oo to t = cxd along the 
upper branch and then returns tot = — oo along the lower branch. Here ipHirt) 
and ■0j(r'i') are the fermion field operators in the Heisenberg picture. 

Parallel to the original procedure of CEID [16] , our main aim is to derive the 
kinetic equation for the one-electron density matrix. Then the key quantity of 
interest is the lesser Green's function G^{rt,r't') = —{ih)~^ (V'jif ('''^')V'-H'('^O) 
since its equal-time value gives the one-electron density matrix: 

Peir,t\r',t) - -ihG<{rt,r't) = (Vj(r'i)V'H(rt)) (3) 

We first derive the equation of motion for the contour-ordered Green's function 
G{rt,r't'). Differentiating G{rt,r't') with respect to time arguments and then 
using equations of motion for the Heisenberg operators, one obtains 

{thdt - h,{r)) Girt, r't') - S{r - r')Sc{t - t') ~ F[r)T^{rt, r't') (4) 



_ h,{r')) Girt, r't') = S{r - r')5c{t - t') - T'^{rt, r't')F{r') (5) 

where h^ir) = +V{r) and 5c{t — t') is the contour delta function [6]. 

Two new nonequilibrium Green's functions are introduced here 

T^{rt,r't') = {ih)-^ {TcXHmH{rt)i;+{r't')) (6) 

T'^{rt,r't') - (zft)-i {Tc^H{rt)i;+,{r't')XH{t')) (7) 
Combining equations (0]) and ^ gives 

Pe = ^[he,Pe]~^[F,tA (8) 

in in 

where we used pe{r,t\r' ,t) ~ —ih\iu\t^t' [{dt + df) G{rt,r't')]^ , and the first 
moment fi{r,t\r',t) is defined as: 

pir,t\r',t) = -ihT<{rt,r't)^-ihT'<{rt,r't) 

= {^^+i{r't)iPH{rt)XH{t)) (9) 

Note that the kinetic equation ([8]) always ensures the conservation of electron 
number, since the relation {N^) — irpe — holds due to the cyclic invariance of 
the trace. 
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We proceed to find the equations of motion for T ^^{rt,r't'): 



{ihdt - h,{r))T^{rt,r't') = 5{r-r')5c{t-t') {XH{t))-F{r)V^,{rt,r't') +—Tx{rt,r't') 

iVl 

(10) 

(-ihdr - K{r'))T^{rt,r't') = S{r - r')6c{t - t') (X^W) - r^,, (rt, r'i')i^(r') 

(11) 

where three new nonequiUbrium Green's functions are defined as: 

Vx{rt,r't') = {ihy^ {TcPH{t)i^H{rt)ij+{r't')) (12) 

T^.irtyt') = (th)-' {TcXl{t)^H{rmti{r't')) (13) 

r;,, {rt, r't') ^ {ih)-^ {TcXHimH{rmUr't')XH{t')) (14) 

In order to obtain a closed set of equations of motion, we decouple the higher- 
order Green's functions T ^^{rt,r't') and T'^^{rt,r't') as follows: 

r^, [rt, r't') « CRR{t)G{rt, r't') (15) 

{rt, r't') w ihD{t, t')G{rt, r't') (16) 

where D{t,t') = {ih)~^ {TcXH{t)XH{t')) is the dressed phonon Green's func- 
tion and CRR{t) = ihD'^{t,t). Using these decoupHng approximations and the 
relation /i(r, t\r' , t) — —ih\\mt->t' [{dt + dt')T^{rt, r't')]^, we combine equations 
lfTO|) and m]) to yield 

where we have applied the Langreth theorem [H [6] to calculate {rt, r't') = 
iHD'^{t,t')G'^{rt,r't'), and the first moment X{r,t\r',t) is defined as 

X{r,t\r',t) = -ihT<{rt,r't) 

= {i>+ir't)Mrt)PH{t)) (18) 

We continue to derive the equation of motion for rx{rt,r't'): 

{ihdt - he{r)) Tx{rt, r't') = S{r - r')Sc{t - t') {Pnit)) - ihKT^{rt, r't') 

+A(ri, r't') - F{r)Tx^{rt, r't') (19) 



{-ihdf ~ K{r')) T^irt, r't') = 5{r - r')5c{t - t') {Pnit)) - T'^^irt, r't')F{r') 

(20) 
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where three new nonequilibrium Green's functions are defined as: 

A(rt,r'i')= J droF{ro){Tc^P+{rot)ibH{rot)i^H{rt)i^jiir't')) (21) 
Tx^{rt,r't') = (ih)-^ {TcPH{t)XHmH{rt)^ti{^'t')) (22) 

T'^^{rt,r't') = izh)-' {TcPHmH{rm+{r't')XH{t')) (23) 

To decouple the above higher-order Green's functions, we make the following 
approximations : 



A(rt,r't')« J droF{ro){^b+irot)^JH{rot)){Tc^PH{rt)^P+{r't')) 
- J dro {ipUrot)ijHirt)) F(ro) {Tci^H{rom+{r't')) 

= ih I droF{ro)pe{ro,t)G{rt,r't')-ih [ droPe{r,t\ro,t)F{ro)G{rot,r't') (24) 



r^^{rt,r't') « {ihy^ {PH{t)XH{t)) (TctPHirmUr't')) 

= CpR{t)G{rt,r't')-'-^Girt,r't') (25) 

r^^(rt, r't') « {TcPH{t)XH{t')) G{rt, r't') (26) 

where CpR{t) = \ {PH{t)XH{'t) + XH{t)PH{t)) ■ Using the above decoupling 
approximations and the relation A(r, t|r', i) — —ihlimt^t' [{dt + dt')T\{rt,r't')]'^ ^ 
we combine equations (fT9|) and (|20|) to arrive at 



A = 4 [he. A] + tr(pei^)Pe + ^ {-F, Pe} - Pe-Fpe - ^CpR [F, p^] - Kp (27) 

in, 2 in. 

where the Langreth theorem is appHed to calculate T'^^{rt, r't') = {XH{t')PH{t)) G'^[rt, r't'). 

Kinetic equations lfT7|) and (|27|) are identical to those derived in Refs. |2l 
I16j . except that an extra term ti{peF)pe appears and a pFp, term disappears in 
equation l(27|) compared with the corresponding equation in the original CEID. 
The reason for the presence of the extra iv{ppF)pe term is that, in the original 
CEID [2 HE], the expansion was with respect to Ai? = R — R, whereas we, for 
convenience, here use X = R — Rq instead. In the original CEID, the pFp. term 
results from higher-order corrections to the Hartree-Fock decoupling for the two- 
electron density matrix in equation l(24|) . which, however is not considered in the 
present formulation. Apart from the two differences, the present formulation of 
CEID is parallel to the original CEID. 



6 



In order to have a closed set of equations, one can derive the perturbation 
expansion for CiiR{t) and CpR{t) by using the expression [4j 



It is obvious that Wick's theorem can be appHed directly to the Tc-products 
of V'-ffo's s-nd V'i^o's- For simplicity, we shall only retain the zero-order term 
in the expansion. Thus we set CiiR{t) = ihD^{t,t) = C (C is a constant) 
and CpR{t) = i {PH„{t)XHoit) + XHo{t)PHo{t)) = 0. We shall hereafter con- 
sider the bare phonon Green's function Do{t, t') only, instead of the dressed one 



The decoupling approximations (fT5|) , lfT6|) , l|24p , (|25l) and (|26l) are the defin- 
ing approximations in the original CEID method [2J. Their key effect is to 
yield single-time equations of motion. These approximations can be well under- 
stood in the framework of diagrammatic perturbation theory as follows. Each 
of them represents a subset of diagrams in the diagrammatic perturbation ex- 
pansion of the corresponding Green's function. As shown in figure 1 where we 
use T'^^ (rt, r't') as an example, the decoupling approximation (fTBl) includes the 
first diagram in figure 1(a) and consequently coincides exactly with the exact 
perturbation expansion at the lowest-order, while it includes only the second 
diagram in figure 1(a) at the second-order. 

In the NEGF-based formulation, CEID can be systematically extended in 
two possible ways. Following the standard equation-of-motion technique, one 
may extend the hierarchy of coupled equations of motion for Green's functions 
and then truncate the hierarchy somewhere by some sort of decoupling proce- 
dure in which a higher-order Green's functions is expressed approximately as 
a product of lower-order Green's functions. Another possibility is to improve 
the decoupling approximations for Green's functions , T'^^ , A, r;^^ and F';^^ 
by adding correction terms. For instance, the last diagram in figure 1(a) (a 
second-order noncrossing diagram in the exact perturbation expansion for F^^ ) 
which is absent from figure 1(b) serves as a second-order correction term to the 
decoupling approximation (flGl) illustrated in figure 1(b). We shall focus on the 
way of making corrections to CEID decoupling approximations, since it does 
not result in a higher hierarchy of coupled equations (see section 13. 2p . 

We have thus rederived the CEID equations of motion from nonequilibrium 
Green's functions and generalized them to an ensemble which allows for a vari- 
able total number of electrons. Moreover, the present formulation allows the key 
approximations in CEID to be quantified in diagrammatic terms and provides 
an in-principle way to improve them. 

3 Comparison with the self-consistent Born ap- 
proximation 

3.1 Weak electron-phonon coupling limit 

The CEID equations of motion to lowest-order in F read: 



{AH{t)BH{t)) (Tc \e-^^c"H,{r). 



D{t,t'). 
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1 r 



he, 11 



(1) 



1 r 
ih 



ih 



c 



M 



(28) 
(29) 



'''' - -h 
ih 



he, A(i)] + tr(p(°)F)p(°) + i {f, } - - if^(i) 



(30) 



where the superscript denotes the order in the couphng strength F. Note that 
the decoupHng approximations (flSl) , (fTBl) , l(24|) , l(25|) and l(26|) are exact to lowest- 
order in i^. Hence the above kinetic equations yield the exact pe . According 
to equation ([3]), the density matrix p'^^ must correspond to the sum of all 
the second-order terms (the Hartree and Fock diagrams) in the perturbation 
expansion of Green's function G(l, 1'), which reads 

G(2)(l,l') = tr(p(")^^) [ d2dt3Do{t2,h)Go{l,2)F{r2)Go{2,l') 



(31) 



Here we use a common short notation (fc) = (rktk). Starting from this equation, 
one can also easily derive equations ((28|) . (|29| and (|30|) with the equation-of- 
motion technique, and furthermore identify that the second term in the right- 
hand side of equation (|30|) comes from the contribution of the Hartree term. 



ih / d2d3Go{l,2)F{r2)Do{t2,t3)Go{2,3)F{r3)Go{3,l') 



Recall that in SCBA the sum of the second-order terms Gg^^j^, i.e. the first 
Born approximation, involves both the Hartree and Fock diagrams and coincides 
exactly with equation (pT]) . Hence, in the weak electron-phonon coupling limit 
CEID agrees exactly with SCBA for an arbitrary nonequilibrium state of the 
electron-ion system, which extends the range of validity of the conclusion in the 
previous comparison [16] for a steady state in the energy domain. 

3.2 The fourth-order in the coupling strength 

We shall go beyond the weak electron-phonon coupling limit and compare CEID 
and SCBA at the fourth-order in F. For simpHcity, Hartree-like diagrams (with 
a closed fermion loop) will be excluded from our analysis. Let us consider the 
fourth-order SCBA Green's function G'^'^^l, 1') = Ai(l, 1') + ^2(1, 1') which is 
represented by two diagrams: 

Ai{l, 1') ^ih f d2d3Go(l, 2)Fir2)Do{h, h)Go{2, 3)F{r3)G'^^^ (3, l') 



^2(1,1')=*^ J d2d3G„{l,2)F{r2)D„{t2,t3)G^^\2,3)F{r3)Go{3,l') 



8 



These diagrams are shown in figure 2(a) and (b) respectively. Here G'^^^(l, 1') = 
ih J d2d3Go(l, 2)F{r2)Da(t2, h)Go{2, 3)F(r3)Go(3, 1') since the Hartree diagram 
is ignored. Kinetic equations for G('*)(1,1') with respect to ti and t'^ yield 

p^e'' = (ifi)^^ he^p^e^ — [iK)^^ [F, /i^'^^j where the corresponding r^'^^(l,l') = 

1') + -82(1, 1') contains two terms: 

B^{l,l')^-ih J d2Do{h,t2)Go{l,2)Fir2)G^^\2,l') 



B2(l,l') = / d2Do{ti,t2)G'^^Hl,2)F{r2)Go{2,l') 



Note that Bi{l, 1') and -82(1, 1') result from the temporal derivatives of Ai(l, 1') 
and ^2(1, 1') respectively. The kinetic equations for F/i "^(l, 1') read 



z;ift,r(f)(l,l') = /ie(ri)r(f)(l,l')-F(ri) [gg(2)(i,i') + 5i(i,i') 



-ih 



M 
(32) 



-zhdt'T^^^l, 1') - r|f)(l, l')h,.{r[) - [ihDoih,t[)G^'\l, 1') + 52(1, l')] Fir[) 

(33) 

where 5i(l,l') = {ih)^ J d2d3i?o(ti, t3)^o(ti, t2)Go(l, 2)F(r2)Go(2, 3)F(r3)Go(3, 1'), 
52(1,1') - (ift)2/d2d3-Do(ti,t2)Go(l,2)F(r2)Go(2,3)F(r3)i?o(t3,ti)Go(3,l'), 
and Tf\l, 1') = Gi(l, 1') + G2(l, 1') where 

Gi(l,l') = -*ft / d2Go(l,2)do(ii,<2)^^(r2)G(2)(2,l') 



G2(l,l') = -*ft / d2do(ti,i2)G(2)(l,2)F(r2)Go(2,l') 



with (io(i,t') = (ih) ^ {TcPHo{t)XHa{t')) ■ The two terms in the square bracket 
in equation (|32| (equation l|33p ') correspond to all the second-order noncrossing 
diagrams in T^, (1, 1') (T'^^ (1, 1')) (see equations JH]) and lO), while GG^^) (1, 1') 
{ihDo{ti,t'i)G^^'' {1,1')) corresponds to the single second-order diagram in the 
decoupHng approximation for (1, 1') (FJ^^ (1' 1')) (see equations (fT5| and lfT6|) ) 
where 5i(l, 1') (^2(1, 1')) is not present. This was illustrated in figure 1, where 
ihDo{ti,t[)G'-^y{l,l') and 52(1,1') in equation (031) are associated with the sec- 
ond and third diagrams in figure 1(a) respectively. Equations ((32| and l(33l) lead 
to 



ih 



he,p' 



(3) 



ih 



C 



f,p: 



(2) 



A(3) 



where Uf,{r,t\r' ,t) = F{r)Sf {rt,r't) - Sf{rt,r't)F{r'). It is seen that H^, 
contributed by the diagrams 5i(l, 1') and 52(1, 1'), serves as a correction to the 



9 



CEID equation of motion for fj,^^^ (cf . equation (fTT)) at the third-order) , so as 
to make CEID equivalent to SCBA at the fourth-order in the coupHng strength. 
However this correction no longer involves just single-time quantities. 

One can proceed to analyze the kinetic equations for (1, 1') in a similar 
manner: 

-ih J drop<-°Hn,h\ro,h)F{ro)G'^^Hroh,r[t[) 

-ih I dropi2)(ri,ti|ro,ti)i^(ro)GoMi,rli;)-zni^r(f)(l,l')(34) 



- ihdt^vfil, 1') = rf (1, l')K{r[) - [zMo(ti, ti)G(2Hl, 1') + ^4(1, 1')] nr'i) 

(35) 

where 53(1,1') = {ihf j d2d3do(ti, i3)^o(ti, ^2)60(1, 2)F(r2)Go(2, 3)F(r3)G'o(3, 1') 
and 54(1, 1') = {ihf J d2dSdo{ti,t2)Goil, 2)F(r2)Go(2, 3)F(r3)i?o(t3, i'i)Go(3, 1'). 
The two terms in the square bracket in equation (|34l) (equation (|35|) ) correspond 
to all the second-order noncrossing diagrams in rA^(l, 1') (r';^^(l, 1')) (see equa- 
tions (I22I) and dSSD), while (Pff„(ii)XH„(ti)) G(2)(1, 1') {ihdo{ti,t[)G'^^\l,l')) 
is the single second-order diagram in the decoupling approximation for rA^(l, 1') 
(r';^^(l, 1')) (see equations {25]) and ^) where ^3(1, 1') (5'4(1, 1')) is not present. 
From equations l(34|) and (|35l) . one obtains 

where Il\{r,t\r' ,t) — F{r)Sf{rt,r't) — Sf{rt,r't)F(r') which serves as a cor- 
rection to the CEID equation of motion for A^'^-' (cf. equation l(27j) at the third- 
order). Note that a term similar to the second term in the right-hand side of 
equation ((3Ql) does not appear because of the exclusion of Hartree-Hke diagrams. 

We have thus identified explicitly the corrections to the CEID equations 
of motion for pS^^ and A*^^-* added by SCBA. These correction terms are con- 
tributed by the diagrams which are absent from the decoupling approximation 
but are present in the complete collection of second-order noncrossing diagrams 
for higher-order Green's functions r^,(l, 1'), 1'), rA;.(l, 1') and r'^^(l, 1')- 

After this correction, CEID becomes equivalent to SCBA at the fourth-order. In 
principle, one may extend this relation between CEID and SCBA to any order 
in F . However, the amount of diagrams increases fast with increasing order so 
that it would not be easy to illustrate their relation in higher order case. 

3.3 Large ionic mass limit 

In the limit of infinite ionic mass, the electron density matrix is determined by 
equation ^ and 
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p=^[K,fi]-^C[F,p,] (36) 

Here the constant C corresponds to the equal-time classical phonon Green's 
function because the oscillator with infinite mass is treated classically (see ap- 
pendix) . In Ref. [16] , it was shown that these coupled equations of motion are 
identical to the corresponding kinetic equations for the following effective elastic 
scattering problem. 

Consider non-interacting electrons linearly coupled to a single infinitely heavy 
classical degree of freedom X, with a distribution xi^) — 5 ~ VC) + 6{X + \/C) 
The system can be described in terms of one-electron density matrix p{X, t) 
which is governed by ihp{X,t) = [h{X), p{X,t)] with the one-electron Hamilto- 
nian h{X) = ho — FX. Define 



p,{t) = J p{X,t)x{X)dX 
ii{t)= f XpiX,t)x{X)dX 



P2{t) = J X^p{X,t)x{X)dX = Cpe 
Then pe{t) is generated exactly by 

ihpe = [ho,pe]-[F,p] (37) 



ihp^[ho,p]-C[F,p^] (38) 

which are identical to equations ^ and (|36|) [16]. We now proceed to solve 
equations l(37|) and (|38l) in integral form. It is proposed that for t > to this 
solution can be written as 

Pe(t) = P<(t,i), p{t)^^i<{t,t) (39) 

where p^{t,t') and /x<(t, i') are defined as 

p<{t,t') = -{ihf j G+it,to)p,{tQ)G^{tQ,t')x{X)dX (40) 

p<{t,t') = -{ihf j XG+(t,to)pe{to)G^{to,t')x{X)dX (41) 

with {ihdt - h{X))G^{t,t') = 6{t-t'). Furthermore G^it,t') can be expressed 
in an iterative form 

G±(t,t') - G^{t,t')-X J G^itX)F{t")G%{t",t')dt" (42) 
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with {ihdt - ho) G^{t, t') = 5(t - t') and F{t) = Oit - to)F. Repeating the use 
of equation (|42| in equations (|40| and l(4T|) , one finds that 

p<it,t') = (1 + G+ • p< (1 + -Eq- • G-) + G+ • E< . G- (43) 
p<{t,t') = -cJ {p<{t,t'')F{t'')G-{t'\t') + G+{t,t'')F{t'')p<{t'',t')) dt" 

= -C j {p<{t,t")F{t")Go{t",t')+G+{t,t")F{t")p<{t",t')) dt" (44) 

with 

E^{t,t')^GF{t)G^{t,t')F{t') 
j:<{t,t')^CFit)p<{t,t')F{t') 

p<{t,t') = -{infG+{t,to)p,{to)Goito,t') 



G^{t,t') = J G^{t,t')x{X)dX = G±(i,0+ j G^{tM)^oitiM)G^{t2,t')dtxdt2 

(45) 

Above B - C stands for A{t,t') = / B{t,t")C{t" ,t')dt" . This convention 
will be used where appropriate hereafter. 

In the infinite mass limit, the CEID equations of motion ^ and (|36l) have 
precisely the same form as the kinetic equations l(37|) and (|38| for the elastic 
scattering problem. Hence, in view of the solution to the elastic scattering 
problem, we may, by analogy, suggest the following ansatz to equations |[8| and 
l(36l) in the context of NEGF (cf. equations l(44l) and (|45l) ): 



Gceid{1A')^Go{1,1')+C j Go{l,2)[F{r2)Go{2,3)F{r3)]G{3,l')d2d3 (46) 
and 

r^(l,l') = -C J Go{l,2)F{r2)G{2,l')d2 

= -C J G{l,2)F{r2)Go{2,l')d2 (47) 

In analogy with the treatment in section ([2]), one can easily verify that 
Pe{r,t\r' ,t) — —ihG'^^jjj{rt,r't) and p,{r,t\r',t) = —ihT'^{rt,r't) are solutions 
to equations ^ and l(36|) . In the large ionic mass limit, the CEID equations 
of motion are thus exactly solvable based on a correspondence between CEID 
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and the elastic scattering problem. Interestingly, the Dyson equation (|46l) is 
consistent with the Born approximation (BA) and Gqeid contains only one 
term at each order in F. 

To compare with SCBA, we need the mixed quantum-classical perturba- 
tion expansion for Gscba which, fohowing the discussion in the appendix, can 
be obtained by replacing the quantum phonon Green's function by the clas- 
sical phonon Green's function ((50|) in the SCBA Dyson equation Gscba = 
Go -I- Go • FDqGF ■ G (the Hartree-hke diagrams are ignored here). In the in- 
finite ionic mass limit, the classical phonon Green's function ijSOl) is a constant 
G. So the SCBA Dyson equation becomes Gscba — Gq + GGq ■ FGF ■ G. 
This equation differs from the CEID solution l(46|) from the fourth-order term 
onwards. For instance, g'^qq^j^ — 2Gp^j^. In the large ionic mass limit, the 
difference between CEID and SCBA is thus precisely quantified. 

4 Conclusions 

In this paper, we have considered, using the nonequilibrium Green's function 
theory, a system of noninteracting electrons linearly coupled to a quantum os- 
cihator. A set of kinetic equations, which determine the one-electron density 
matrix, are derived with the equation-of-motion technique. Our work estab- 
hshes a rigorous connection between CEID and NEGF, and extends the scope 
of CEID to a general nonequilibrium ensemble that allows for a variable total 
number of electrons. By perturbation theory, the decoupling approximations 
used in the CEID methodology can be quantified in diagrammatic terms. 

We have compared the limiting behavior of CEID and SCBA analytically. 
In the weak electron-phonon coupling limit, they agree exactly for a general 
nonequilibrium state of the system. In the large ionic mass limit, where CEID 
corresponds to an elastic scattering problem and can be solved exactly, the 
difference between CEID and SCBA emerges from the fourth-order term and 
can be quantified. In particular, we ihustrate the connection between CEID 
and SCBA at the fourth-order in the couphng strength. We find that, CEID 
occupies a special place between BA and SCBA, such that CEID is simpler than 
SCBA but is an improvement over BA, in that CEID conserves total number of 
electrons. The lowest-order SCBA corrections to the CEID equations of motion 
(see section [3^ no longer involve just single-time quantities. This illustrates 
the sense in which CEID can be thought of as the simplest particle-number 
conserving approximation that, in addition, retains just single-time quantities. 

The present formulation of CEID can be extended to include multiple quan- 
tum oscillators. The purpose of the single-oscillator model calculation is to 
illustrate the analytical features of CEID in a simple way so that an analytical 
comparison of CEID and SCBA can be made. Like SCBA, the present method 
for CEID is not applicable to problems with strong electron-phonon correlations, 
which, however, have been addressed by another CEID scheme [1^ recently. 
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Appendix 

In this appendix, we shall consider the Green's function for a system of quan- 
tum electrons coupled to a classical oscillator and then show how to develop a 
diagrammatic perturbation expansion for it. The system Hamiltonian takes the 
form id]) but X and P are classical variables now. The mixed quantum-classical 
Green's function, as usual, can be written as 

G(l,l') = m-'(Tc [*^„(l)vI/+^(l')e-^^^^-oM'^-]) 

where the angular bracket (• • •) ~ J dXdPtr {po ■ ■ •)■ To evaluate this Green's 
function, we must provide a procedure for evaluating the average of products 
of classical coordinates AT's, while the average of products of electronic field 
operators can be evaluated by Wick's theorem. 

Consider a classical oscillator with position, momentum and energy 

X{t) = A cos {cut - (j)) 



P(t) = -AloM sin {ujt - (j)) 



+ -Muj^X^ = -Muj^A^ 

2M 2 2 

sampled from the canonical distribution 

2n 

Let (. . .) denote averaging over p. Changing variables from {X, P) to {A, (j)) 
with dXdP MLoAdAdcj), 




Our aim is to establish the relation 
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L = {X{h)X{t2) . . . X{hN)) - i? = {X{ti)X{t2)) {X{t3)X{U)) . . . {X{t2N-l)X{t2N)) 



+ all other pairings (48) 

Write 

X{u) = ^ (e'M--*) + e-'M--*)) = ^^-t^, a, = Ae'M.^0) 

Consider L. Expand and integrate over (f). Only terms with N a's and TV a*'s 
survive. There are 

^ 7V!iV! 

such terms. Each is of the form 

2^ (/3Mc^2)iv ^ e J ^49j 

with one such term occurring in L for each of the Nl possible groupings of the 
2N indices into two groups of N, {{ki . . . fcjv) , (fcAr+i . . . k2N)}- 

Now consider R. Note that the classical phonon Green's function 

{XiU)Xit,)) -^cos.iU^t,) = ^ (e^^(.-.) +e-(-*^)) (50) 

In R there are 

(2NV 

iVp = (27V-l)(2iV-3)...l = i^ 
different pairings. Each pairing contributes 2^ terms, each of the form 



1 I /'p«^(tfci+---+tfcjv) X e"'"(*'"«+i+---+*'=2jv)^ (fi-i) 

2^ (BMiv^W V ^ J ^^^^ 



2^ {/3MLu^y 
Thus, R is composed of 

(2jV)!^^_ {2Ny. 

2^Ni m 



terms, each of the form l(5T|) . 

By symmetry, every grouping of indices {(fci . . . kpf) , (fcjv+i . . . k2N)} that 
occurs in L occurs in R and vice versa. Further, by symmetry, if a given 
grouping {(fci . . . k]^) , (fcAr+i . . . fc2Ar)} occurs G times in R, then every other 
grouping must also occur G times in R. There are Nl groupings. Hence 
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and thus every grouping {(fci . . . fcjv) , {kN+i ■ ■ ■ ^Af)} occurs A^! times in R. 
Using this and (|5T|l . we see that in R each distinct grouping occurs with a 
prefactor 

1 iV! 

2^ (/3Mw2)Af 

which is the same as the prefactor with which each grouping occurs in L (see 
equation l(49|)). Hence L = R. This relation (Wick's theorem) allows us to eval- 
uate the mixed quantum-classical Green's function as a perturbation expansion 
involving only wholly paired nuclear coordinates. 
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Figure captions 



Figure 1. (a) Zero- and second-order noncrossing diagrams in the exact F^^^ (rt, j-'t') 
(see equation mi)), (b) The diagrammatic representation of the CEID 
decoupUng approximation for T'^^{rt,r't') (see equation (flGl) ). The thick 
(thin) straight line represents dressed (bare) electron Green's function. 
The thick (thin) wavy line represents dressed (bare) phonon Green's func- 
tion. We ignore Hartree-like diagrams and corrections to a phonon line. 

Figure 2. The diagrammatic representation of the fourth-order SCBA Green's 
function G('')(l,l'). 
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Y. Wang, Figure 1 
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Y. Wang, Figure 2 
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